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ABSTRACT 

The initial-final mass relationship connects the mass of a white dwarf with the mass 
of its progenitor in the main-sequence. Although this function is of fundamental im- 
portance to several fields in modern astrophysics, it is not well constrained either from 
the theoretical or the observational points of view. In this work we revise the present 
semi-empirical initial-final mass relationship by re-evaluating the available data. The 
distribution obtained from grouping all our results presents a considerable dispersion, 
which is larger than the uncertainties. We have carried out a weighted least-squares 
linear fit of these data and a careful analysis to give some clues on the dependence 
of this relationship on some parameters such as metallicity or rotation. The semi- 
empirical initial-final mass relationship arising from our study covers the range of 
initial masses from 1.0 to 6.5 Mq, including in this way the low-mass domain, poorly 
studied until recently. Finally, we have also performed a test of the initial-final mass 
relationship by studying its effect on the luminosity function and on the mass distribu- 
tion of white dwarfs. This was done by using different initial-final mass relationships 
from the literature, including the expression derived in this work, and comparing the 
results obtained with the observational data from the Palomar Green Survey and the 
Sloan Digital Sky Survey (SDSS). We find that the semi-empirical initial-final mass 
relationship derived here gives results in good agreement with the observational data, 
especially in the case of the white dwarf mass distribution. 

Key words: white dwarfs — stars: evolution, luminosity function, mass function — 
open clusters and associations: general. 



1 INTRODUCTION 

Tlie initial-final mass relationship of white dwarfs links the 
mass of a white dwarf with that of its progenitor in the 
main-sequence. This function is of paramount importance 
for several aspects of modern astrophysiscs such as the de- 
termination of the ages of globular clusters and their dis- 
tances, the study of the chemical evolution of galaxies, and 
also to understand the properties of the Galactic population 
of white dwarfs. However, we still do not have an accurate 
measurement of this relationship and, consequently, more 
efforts are needed from both the theoretical and the obser- 
vati onal perspect i ves to improve it. 

IWeidemannI (Il977l ) carried out the first attempt to 
empirically m ap this relationsh ip, and provided also a re- 
cent revision (|Weidemannll2000l ). Although many improve- 
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ments have been achieved in these 30 years, there are still 
some pieces missing in the puzzle. For instance, the depen- 
dence of this function on different parameters is still not 
clear (e.g. metallicity, magnetic field, angular momentum). 
On the other hand, numerous works have dealt with the 
calculation of a theore t ical initial-fin al mass relationship 
(|Dommguez et al.|[l999l : iMarigol bOOll ). but the differences 
in their evolutionary codes, such as the treatment of con- 
vection, the value of the assumed critical mass, which is the 
maximum mass of a white dwarf progenitor, or the mass loss 
prescriptions used lead to very different results. The main 
differences between the different theoretical approaches to 
the initia l-final mass relation ship have been extensively dis- 
cussed in IWeidemannI (I2OO0I ). 

From an observational perspective, most efforts up to 
now have focused on the observation of white dwarfs in open 
clusters, since this allows to infer the total age and the 
original metallicity of white dwarfs belonging to the clus- 
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ter I Williams et al.ll2004l : iKalirai et aLlbOOSh . Open clusters 
have made possible the derivation of a semi-empirical initial- 
final mass relationship using more than 50 white dwarfs, 
although only covering the initial mass range between 2.5 
and 7.0 Mq because stellar clusters are relatively young 
and, hence, the white dwarf progenitors in these clusters 
are ge nerally massive. The recent study of iKalirai et al.l 
l|200^ based in old open clusters (NGC 7789, NGC 6819 
and NGC6791) has extended this mass range to smaller 
masses. A parallel attempt to cover the low-mass range of 
the initial-fina l mass relationship has been carried out by 
ICatalari" et al.l l|2008l ). This was the first study of this rela- 
tionship b ased in common pro per motion pairs. The stars 
studied in ICatalan et al.l (|2008l ) are at shorter distances in 
comparison with star clusters and this allows a better spec- 
troscopic study of both members of the pair, obtaining their 
stellar parameters with accuracy. At the same time, the 
study of these pairs enables a wide age and metallicity cov- 
erage of the initial-final mass relationship. 

The aim of this work is to perform a revision of the 
initial-final mass relationship taking into account our recent 
results fro m studying white dw arfs in common proper mo- 
tion pairs ijCatalan et al.|[2008l ) and the available data that 
are currently being used to define the initial-final mass rela- 
tionship. The paper is organized as follows. In §2 we present 
the analysis of the data that we will use to define the semi- 
empirical initial-final mass relationship. Section 3 is devoted 
to the detailed analysis of the semi-empirical initial-final 
mass relationship derived in this paper and to give some 
clues on its dependence on some parameters, such as metal- 
licity. In §4 and §5 we compute the luminosity function and 
the mass distribution of white dwarfs considering difi'erent 
initial-final mass relationships and compare our results with 
the available observational data. Finally in §6 we summarize 
our main results and we draw our conclusions. 



2 ANALYSIS OF CURRENT AVAILABLE 
DATA 

2.1 Open clusters and visual binaries 

We have carried out a re-analysis of the available data cur- 
rently used to define the semi-empirical initial-final mass 
relationship, which is mainly based on white dwarfs in open 
clusters. We have used the white dwarf atmospheric param- 
eters (Toff and logjr) derived by other authors, as well as 
the ages and metallicities of the clusters reported in the 
literature. To obtain the final and initial masses w e fol- 
lowed the procedure described in lCatalan et al.l l|2008l ). This 
procedure consists in deriving the final mass (Mf) and the 
cooling time of each white dwarf from the atmos phe ric pa - 
rameters and the cooling sequences of i Salaris et all l|2000l ). 
These cooling tracks consider a carbon-oxygen (CO) core 
white dwarf (with a larger abundance of O at the centre 
of the core) with a H thick envelope ontop of a He bufi'er, 
g(H) = A/h/M = IQ-'' and q(He) = Mho/M = 10"^ As it 
will be shown in the next section, the thicknesses of these 
envelops are very important in the cooling of white dwarfs. 
These improved cooling sequences include an accurate treat- 
ment of the crystallization process of the CO core, including 
phase separation upon crystallization, together with up-to- 



date input physics suitable for computing white dwarf evo- 
lution. Since we know the total ages of these white dwarfs 
(from the age of the cluster) we derived the main-sequence 
lifetimes of the progenitors, and from these, their initial 
masses using the stellar tracks of iDomfnguez et al.l l|l999t ). 
At present, the atmospheric parameters of white dwarfs can 
be determined with accuracy if they correspond to the DA 
type, that is, if their spectra shows uniquely the hydrogen 
absorption lines. For this reason, and in order to keep con- 
sistency in the cooling sequences used we only consider DA 
white dwarfs in our study. In Table [T] we give the initial and 
final masses that we have recalculated for white dwarfs in 
open clusters. Other information such as the atmospheric 
parameters, cooling times, main-sequence lifetimes of the 
progenitors and metallicities are also given. Those white 
dwarfs for which we have obtained a cooling time longer 
than its total age have not been included in our study, since 
this is a good indication that the star does not belong to 
the cluster, or it has not a CO core. To compute the errors 
of the final masses we have taken into account the errors of 
the atmospheric parameters. In the case of the initial masses 
we have taken into account the errors of the cooling times, 
which come from the atmospheric parameters, and the er- 
rors of the total ages of the white dwarfs. When no error 
for the total age was given in the literature, a value of 10 
per cent was adopted. Finally, as reported in lCatalan et al.l 
(2008), to derive the initial masses of white dwarfs belonging 
to common proper motion pairs, the error of the metallicity 
was also taken into account. 

NGC 2099 (M37) 

IKalirai et al.l (|2005l ) performed spectroscopic o bservations of 
30 wh ite dwarfs belonging to NGC 2099 (M37'). IKalirai et all 
l|200ll ') determined the age and metallicity of M37, 650 Myr 
and Z = 0.011, respectively. We assume an error of 10 per 
cent in the age of M37. 

NGC 2168 (M35) 

We use the atmospheric parameters reported by W il liams , 
Bolte & Koester (2004). iBarrado v Navascues et al.l (|200ll ) 
estimated the metallicity of this cluster, [Fe/H]= —0.21 ± 
0.10. Thus, we use the stel l ar tra cks corresponding to Z — 
0.012. As in Wi lliams et al.l l|2004[ ) we use the age derived by 
,von Hippcl (20051 ). 150 ± 60 Myr. The atmospheric param- 
eters oMjAWDSland_JjAWDS27 correspond to new data 
from lFerrario et all (120051 ). 

NGC 3532 

The cluster data are from iKoester fc Reimerd (Il993|). al- 
thoug h we have used the latest resul ts of iKoester fc ReimersI 
(1 19961) report e d in iFerrario et al.l (|2005l ) . According to 
iFerrario et ai] (|2005l ) the age of this cluster is 300±150, 
and the metallicity is [Fe/H]= —0.022 (Twarog, Ashman 
fc Anthony- Twarog, 1997) or [Fe/H]= -0.02 (Chen, Hon fc 
Hang, 2003), so we have used the tracks corresponding to 
Z = 0.019. We ha ve also included the data for three white 
dwarfs reported in iReimers fc Koesteil (|l989D . although the 
resolution is a bit lower. 
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Table 1. Rosiilts from rc-cmihiatiiig tho a,\'aila})lc data. 



WD 


Teff (K) 


lo£ fdexl 


Mr (Mn,) 


t^-nni fGvr) 

'^COOl / 


trtrnv fGvr) 


Mi (Mn.) 


Z 


NGC 2099 (M37) 














0.011 


WD2 


19900±900 


8.11±0.16 


0.69±0.07 


0.093±0.019 


0.56±0.07 


r, 7r,-|-0.12 

^•'^-0.10 




WD3 


18300±900 


8.23±0.21 


0.76±0.09 


0.152±0.034 


0.50±0.07 


r. 00-1-0. 16 

^•8^-0.13 




WD4 


16900±1100 


8.40±0.26 


0.86±0.12 


0.259±0.072 


0.39±0.10 


3.06+SSo 




WD5 


18300±1000 


8.33±0.22 


0.82±0.11 


0.192±0.078 


0.46±0.10 


2.90+°?^ 

— 0.19 




WD7 


17800±1400 


8.42±0.32 


0.88±0.14 


0.219±0.082 


0.43±0.10 


2.96+°-30 

— \J.21 




WD9 


15300±400 


8.00±0.08 


0.61±0.03 


0.182±0.015 


0.47±0.07 


r, 00+0. 16 

— 0.13 




WDIO 


19300±400 


8.20±0.07 


0.74±0.03 


0.120±0.010 


0.53±0.06 


2.76+°" 

— U. 11 




WDll 


23000±600 


8.54±0.10 


0.98±0.04 


0.136±0.014 


0.51±0.07 


2.80+0" 

— 0.11 




WD12 


13300±1000 


7.91±0.12 


0.56±0.05 


0.239±0.043 


0.41±0.08 


3.01+°?2 

— 0.17 




WD 13 


18200±400 


8.27±0.08 


0.78±0.03 


0.167±0.016 


0.48±0.07 


2.86«;}^ 




WDM 


11400±200 


7.73±0.16 


0.45±0.07 


0.282±0.039 


0.37±0.07 


1 9+0.26 
■^•-^^-0.19 




WD16 


13100±500 


8.34±0.10 


0.82±0.05 


0.480±0.062 


0.17±0.09 


4.10+llt 
^■^^—0.57 




NGC 2168 (M35) 














0.012 


LAWDS 1 


32400±512 


8.40±0.12 


0.89±0.06 


0.023±0.006 


0.127±0.060 


4.61+?,1! 




LAWDS 2 


32700±603 


8.34±0.08 


0.85±0.04 


0.017±0.003 


0.133±0.060 


~l'21 
— 0.60 




LAWDS 5 


52600±1160 


8.24±0.09 


0.82±0.04 


0.0022±0.0001 


0.148±0.060 


4.35+°-?! 

— 0.52 




LAWDS 6 


55200±897 


8.28±0.06 


0.84±0.03 


0.0020±0.0001 


0.148±0.060 


4 35+0.98 




LAWDS 15 


29900±318 


8.48±0.06 


0.94±0.03 


0.046±0.005 


0.104±0.060 


4 qq+2-25 




LAWDS 22 


54400±1203 


8.04±0.12 


0.72±0.04 


0.0025±0.0002 


0.147±0.060 


4.35+°-?^ 




LAWDS 27 


30500±397 


8.52±0.06 


0.98±0.03 


0.048±0.004 


0.102±0.060 


2 37 
5.03 n QQ 




NGC 3532 














0.019 


3532- WD 1 


28000±2000 


8.45±0.45 


0.92±0.20 


0.055±0.036 


0.245±0.154 


3.70+,V?5 

— 0.55 




3532- WD5 


28500±2000 


7.8±0.3 


0.55±0.11 


0.013±0.002 


0.287±0.150 


KO+1-05 
— 0.4d 




3532- WD6 


28500±3000 


8.5±0.5 


0.96±0.23 


0.060±0.046 


0.240±0.157 


70+1. 80 
— 0.57 




3532-WD8 


23367±1065 


7.71±0.15 


0.48±0.05 


0.023±0.003 


0.277±0.150 


— 0.48 




3532- WD9 


29800±616 


7.83±0.23 


0.56±0.08 


0.012±0.001 


0.288±0.150 


51 -1-1.03 

— 0.46 




3532- WDIO 


19267±974 


8.14±0.27 


0.71±0.11 


0.112±0.030 


0.188±0.153 


4.05+^-^? 

— 0.73 




Praesepe 














0.027 


WD0836+201 


16629±350 


8.01±0.05 


0.62±0.02 


0.144±0.009 


0.481±0.051 


2 97+011 

^•»'_0.10 




WD0836-I-199 


14060±630 


8.34±0.06 


0.82±0.03 


0.39±0.04 


0.235±0.064 


3.76+°1E 

^ ' ^ — 0.28 




WD0837-M99 


17098±350 


8.32±0.05 


0.81±0.03 


0.22±0.02 


0.405±0.054 


ic:+0.15 
— 0.13 




WD0840-f-200 


14178±350 


8.23±0.05 


0.75±0.02 


0.31d=0.02 


0.315±0.054 


3.42+"-2l 

— (J. 17 




WD0836-I-197 


21949±350 


8.45±0.05 


0.91±0.03 


0.13±0.01 


0.495±0.051 


2 94+"- 11 
— 0.09 




WD0837-I-185 


14748±400 


8.24±0.055 


0.76±0.02 


0.288±0.018 


0.337±0.053 


0C-+O.2O 
— 0. lb 




WD0837-I-218 


16833±254 


8.39±0.03 


0.85±0.01 


0.251±0.009 


0.374±0.051 


90+0. 16 

* ZO C\ TO 

— 0.13 




WD0833-I-194 


14999±233 


8.18±0.035 


0.72±0.01 


0.246±0.009 


0.379±0.051 


3.22+0" 

— 0.13 




WD0840+190 


14765d=270 


8.21±0.03 


n.74±0.01 


0.27±0.01 


0.355±0.051 


3 29+0-17 




WDO8IO+2O0 


14527±282 


8.24±0.04 


0.76±0.02 


n.3n±o.oio 


0.325±0.052 


Q on + 0.2() 

-O.Ki 




\\l)()si:i-l8l 


llil),s±202 


,S.22±0.0i 


().7r)±().n2 


().2l)r)±0,012 


().:i:-!0±O.Or)l 


0,1:, 




Hyades 














0.027 


WD0352+098 


14770±350 


8.16±0.05 


0.71±0.02 


0.25±0.01 


0.375±0.051 


o 90-1-0. 16 

■^•^■^-o.is 




WD0406-I-169 


15180±350 


8.30±0.05 


0.79±0.02 


0.29±0.02 


0.335±0.054 


cic:+0.20 
■^•■3O_0.16 




WD0421-I-162 


19570±350 


8.09±0.05 


0.68±0.02 


0.096±0.006 


0.529±0.050 


9 8S+010 




WD0425-I-168 


24420±350 


8.11±0.05 


0.70±0.02 


0.038±0.003 


0.587±0.050 


9 70-1-0.09 
^- '"-0.08 




WD0431-I-125 


21340±350 


8.04±0.05 


0.65±0.02 


0.060±0.004 


0.565±0.050 


9 09+0.09 
^•"^-0.08 




WD0438+108 


27390±350 


8.07±0.05 


0.68±0.02 


0.018±0.001 


0.607±0.050 


9 7c:+0.08 




WD0437-I-138 


15335±350 


8.26±0.05 


0.77±0.02 


0.26±0.01 


0.365±0.051 


9fi+017 
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Table 1. Results from re-evaluating the available data (continued). White dwarfs in common proper motion pairs are listed in the last 
section of this table (CPMPs) 



WD 


T,s (K) 


logg (dex) 


Mi (Mq) 


icooi (Gyr) 


tprog (Gyr) 


Mi (Mq) 


z 


NGC 2516 














0.02 


2516- WDl 


28170±310 


8.48±0.17 


0.95±0.08 


0.060±0.012 


0.081±0.012 


= =4+0.39 
^•°^-0.29 




2516- WD2 


34200±610 


8.60±0.11 


1.01±0.04 


0.035±0.008 


0.106±0.008 


5 nr,+0.14, 
''•'-'•^-0.13 




2516- WD3 


26870±330 


8.55±0.07 


0.99±0.03 


0.082±0.005 


0.059±0.005 


6.44+0-32 

— 0.29 




2516- WD5 


30760±420 


8.70±0.12 


1.07±0.05 


0.074±0.014 


0.067±0.014 






NGC 6791 














0.040 


WD7 


14800±300 


7.91±0.06 


0.56±0.02 


0.17±0.01 


8.33±0.85 


1.086+!',"« 

—0.038 




WD8 


18200±300 


7.73±0.06 


0.47±0.02 


0.063±0.005 


8.437±0.850 


1 nsi +0.044 




NGC 7789 














0.014 


WD5 


31200±200 


7.90±0.05 


0.60±0.02 


0.010±0.0002 


1.39±0.14 


1 84+°-"'^ 




WD8 








n 029+0 (193 


1.37±0.14 


1 

^•"''-0.05 




WD9 


20900±700 


7.84±0.12 


0.54±0.04 


0.042±0.006 


1.36±0.14 


1 8fi+0.08 
1-^0-0.05 




NGC 6819 














0.017 


WD6 


21100±300 


7.83±0.04 


0.54±0.02 


0.041±0.002 


2.46±0.25 


-^•'"-0.05 




WD 7 


16000±200 


7.91±0.04 


0.57±0.01 


0.139±0.005 


2.36±0.25 






Pleiades 














0.019 


WD0349+247 


32841±172 


8.63±0.04 


1.03±0.02 


0.048±0.004 


0.071±0.008 


C 07 + 0.31 




Sirius B 


25000±200 


8.60±0.04 


l.OitO.Ol 


0.108±0.003 


0.129±0.012 


^•°'-0.16 


0.020 


CPMPs 
















WD0315-011 
WD0413-077 


7520±260 
16570±350 


8.01±0.45 
7.86±0.05 


0.60±0.20 
0.54±0.02 


1.20±0.56 
0.112±0.008 


2.97^3;?^ 
0.96±0.37 


-^•^"-0.28 

r) rir7 + 0.53 

Z.U(_Q 27 


0.016+0.003 
0.008+0.001 


WD1354+340 
WD1544-377 


13650±420 
10600±250 


7.80±0.15 
8.29±0.05 


0.50±0.04 
0.78±0.02 


0.20±0.02 
D.76±D.05 


0.18±0.50 


1 46+°-3l 
-^•^'-'-o.oo 

4 13+- 
^•^•-'-1.49 
3 45+0.65 
■^•^'^-0.35 


0.015+0.002 
0.021+0.003 


WD1620-391 


24900±130 


7.99±0.03 


0.63±0.01 


0.026±0.001 


0.30±0.12 


0.020+0.003 


WD1659-531 


14510±250 


8.08±0.03 


0.66±0.01 


0.24±0.01 


r, 97+0.34 
-0.32 


1 c:o+0.08 

-■-•^S-o.os 


0.019+0.004 



Hyades and Praesepe 

We have used the data of lClaver et all (|200ll ) for the Hyades 
and part of the Pra esepe sample. We h ave included also 
two stars studied by iDobbie et al.l (l2004l) an d some recent 
results on s i x new stars ( Dobbie et al.ll2006l) . According to 
Ivon Hippel (120051 ) . the Hyades cluster has an age of 625±50 
Myr. We assume the same value for Praesepe, since both 
belong to the same Hyades supercluste r . The metallicity is 
[Fe/H]= 0.13 according to IChen et"al] l|2003l ). so we used 
the stellar tracks corresponding to Z — 0.027. 



NGC 2516 

We consider the atmo spheric parameters derived by 
iKoester fc Reimersl(|l996l ). The age of the cluster is that de- 
termined by Meynet, Mermilliod & Maeder (1993), 141±2 
Myr. This cluster has a solar metallicity according to Jef- 
fries, James & Thurston (1998). 



NGC 6791, NGC 7789 and NGC 6819 

We in clude the recent results obtained by iKalirai et alJ 
(i2008h based on spectroscopic observations of white dwarfs 
in different old clusters. These re sults, as well as the ones 
obtained by ICatalan et al are somewhat relevant, 

since constitute the first constraints on the low-mass end 
of the initial-final mass relationship. NGC 6791 is one of 
the oldest and mos t metal-rich open clusters. According to 
iKalirai et all (|2007l ). NGC 6791 has an age of 8.5 Gyr and 
[Fe/H]= -1-0.3 ± 0.5. We have used the stellar tracks cor- 
responding to Z = 0.04. The white dwarfs of this cluster 
have relatively low masses (< 0.47 Mq), so the majority 
of white dwarfs belonging to this cluster could have a He 
core (Kalirai ot al. 2007). Although it has been suggested 
that He-core white dwarfs could be formed also by single- 
evolution (Kilic, Stanek & Pinsonneault, 2007), we prefer 
not to take these white dwarfs into account, and only con- 
sider the two confirmed cluster members with masses above 
this value. In this case, we assume an error of 10 per cent in 
the age of NGC 6791. 

The white dwarf population of NGC 7789 and NGC 
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6819 were also studied recently by iKalirai et al.l(|2008l ). Ac- 
cording to their results, NGC 7789 has an age of 1.4±0.14 
Gyr and a metaUicity oi Z = 0.014, while NGC 6819 has an 
age of 2.5±0.25 Gyr and Z = 0.017. 



Pleiades 

WD0349-I-247 is the only known Pleiades white dwarf. 
We have us e d the atmospheric parameters derived by 
iDobbie et all l|2006l ). The age of the Pleiades is 125±8 Myr 
according to Rebolo, Martin & Magazzu (1992) a nd it has 
a metaUicity of [Fe/H]= -0.03 (jChen et al.ll2003h . We use 
the stellar tracks corresponding to Z — 0.019. 



Sirius B 

Sirius B is a very well-known DA white dwarf that belongs 
to a visual binary system. Since the separation between the 
members of this system is ~ 7 AU it can be assumed that 
there has not been any significant interaction between them 
during th e AGB phase of the progenitor of the white dwarf 
member (|Liebert et al.|[2005ah . Thus, it can be considered 
that Sirius B has evolved as a single star, and therefore it 
is a good candidate for studying the initial-final mass rela- 
tionship. In t his work, we u s e the atmospheric parameters 
determined bv lLiebert et al.1 l|2005ah . an age of 237±12 Myr 
and solar metaUicity. 



2.2 Globular clusters 

In principle, globular clusters could also be used for improv- 
ing the initial-final mass relationship since they have some 
characteristics which are similar to those of open clusters, 
but we have decided not to include them since the study of 
white dwarfs in globular clusters still suffers from large un- 
certainties. Globular clusters contain thousands to millions 
of stars. However, white dwarfs belonging to them are usu- 
ally very faint due to their large distances, which difficults 
obtaining high signal-to-noise spectra. Moreover, since glob- 
ular clusters usually have crowded fields it is also difficult to 
isolate each star properly when performing spectroscopic ob- 
servations. Therefore, the reduction procedure, mainly the 
background substraction, is also more complicated. Up to 
now, the white dwarf population of globular clusters has 
been mainly used to determine their distances and ages 
l| Richer et al .11200^ . An important attempt to derive the 
masses of white dwarfs in globular clusters (NGC 6397 and 
NGC 6752) was carried out by Mochlcr ct al. (2004). How- 
ever, due to the low signal-to-noise and resolution of the 
spectra, they were not able to derive their spectroscopic 
masses independently and had to use a combination of spec- 
tra with photometry. Although they calculated an average 
mass for the white dwarfs in NGC 6752 we have decided not 
to include this value in our work. Instead, we prefer to use 
data of individual white dwarfs rather than a bin ned value 
for a cluster, contrary to that done in other works l|Williamsl 
l2007l : lKalirai et al.ll2008l ). 
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Figure 1. Final masses versus initial masses of the available clus- 
ter and common proper motion pairs data. 



2.3 Common proper motion pairs 

In the case of the common proper motion pairs, the proce- 
dure that we followed to derive the final and initial masses 
of the white dwarfs is explained in detail in iCatalan et al.l 
(|2008h . It mainly consisted in performing independent spec- 
troscopic observations of the components of several common 
proper motion pairs composed of a DA white dwarf and a 
FGK star. From the fit of the white dwarf spectra to syn- 
thetic models we derived their atmospheric parameters, and 
from these t heir masses and coo ling times using the cooling 
sequences of ISalaris et al] (|2000l ). In order to derive the ini- 
tial masses of the white dwarfs we performed independent 
high-resolution spectroscopic observations of their compan- 
ions (FGK stars). Since it can be assumed that the mem- 
bers of a common proper motion pair were born simulta- 
neously a nd with the s ame chemical composition (|Wegnerl 
1973; Os walt et al.l |l988), the ages and metallicities of both 
members of the pair should be the same. From a detailed 
analysis of the spectra of the companions we derived their 
metallicites. Then, we obtained their ages using either stel- 
lar isochrones, if the star was moderatly evolved, or the 
X-ray luminosity if the star was very close to the ZAMS 
(Ribas et al. in preparation). Once we had the total age 
of the white dwarfs and the metaUicity of their progenitors 
we derived their initial masses using the stellar tracks of 
'Dommgu ez et al.1 ([l99i). In Table [1] we give the initial and 
final masses resulting from this study. 



3 THE INITIAL-FINAL MASS RELATIONSHIP 

In Fig. [1] we present the final versus the initial masses ob- 
tained for white dwarfs in common proper motion pairs and 
open clusters. The observational data that can be used to 
define the semi-empirical initial-final mass relationship con- 
tains now 62 white dwarfs. It is important to emphasize that 
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all the values below 2.5 Mq correspond to our data obtained 
from common pr oper motion pairs (C PMPs) and the recent 
data obtained bv lKalirai et alj ()2008l ) — K08. Before these 
studies no data for these small masses were available, since 
white dwarfs in stellar clusters are usually more massive, 
especially if the clusters are young. The coverage of the low- 
mass end of the initial-final mass relationship is specially 
important since it guarantees, according to the theory of 
stellar evolution, the study of white dwarfs with masses near 
the typical values, M ~ 0.57 Mq, whic h represent about 9 
per cent of the white dwarf population l|Kepler et al.|[2007l ). 
Thus, these new data increase the statistical significance of 
the semi-empirical initial-final mass relationship. 

A first inspection of Fig. [1] reveals that there is a clear 
dependence of the white dwarf masses on the masses of their 
progenitors. In Fig. [T] we have als o plotted the the o retica l 
initial-final mass relationships of Dom mguez et all (Il999l ) 
for difi^erent metallicities to be consistent with the stellar 
tracks used to derive the initial masses. Although the distri- 
bution presents a large dispersion, a comparison of the ob- 
servational data with these theoretical relationships shows 
that they share the same trend. However, it should be noted 
that for each cluster the data presents an intrinsic spread in 
mass. The dispersion varies from cluster to cluster, but it is 
particularly noticeable for the case of M37. Nevertheless, it 
should be taken into account as well that the observations 
of M37 were of poo rer quality than the rest of the data 
l|Ferrario et al.ll2005h . 

3.1 Main systematic uncertainties 

The results obtained are dependent on different assumptions 
and approaches that we have considered during the proce- 
dure followed to derive the final and initial masses, we dis- 
cuss them separately. 

3.1.1 Thicknesses of the H and He envelopes 

The fact that the observed white dwarf masses in clusters 
scatter considerably in the same region of initial m asses, 
as poi nted out by other authors — see, for instance, iReidI 
1 19961 ) — rnay indicate that mass loss could depend more 
on individual stellar properties than on a global mechanism, 
and that it could be a stochastic phenomenon, especially on 
the AGB phase. Mass loss has a large impact on the final 
composition of the outer layers of white dwarfs, since it de- 
fines the thicknesses of the outer He and H (if present) layers. 
In fact, another reason that may explain why white dwarfs 
with different final masses could have progenitors with very 
similar initial masses is the assumption of a given inter- 
nal composition and outer layer stratification of the white 
dwarfs under study. The thicknesses of the H and He lay- 
ers is a key factor in the evolution of white dwarfs, since 
they control the rate at which white dwarfs cool down. In 
this work we have used cooling sequences with fixed thick- 
nesses of these envelopes, which might be more appropriate 
in some cases than in others. In fact, the exact masses that 
the layers of H and He may have is currently a matter of 
debate being the subject o f seve ral studies. For instance, 
iPrada Moroni fc Stranierol (|2002l ) computed models reduc- 
ing the thickness of the H envelope to 5(H) — 2.32 x 10~®, 



obtaining cooling times shorter than those obtained in this 
work assuming a thick envelope (q(H) = 10~^). This is natu- 
ral since H has a larger opacity than He, and H is the major 
insulating component of the star. In the case of an even 
thinner H envelope (g(H) — 10^^") the cooling age could 
be reduced in 10 per cent (1 Gyr) at log(L/LQ) = —5.5. 
Thus, the uncertainty in the cooling times could be rele- 
vant in some cases, which would affect the estimates of the 
progenitor lifetimes and in turn, the initial masses derived. 

In order to estimate the effect that the thicknesses of 
H and He envelopes may have in the initial masses derived 
here, we have repeated the calculations of section 2 but us- 
ing the cooling sequences of Fontaine, Brassard & Bergeron 
(2001) for a 50/50 CO core white dwarf with a standard He 
envelope, g(He) = 10~^, and two different thicknesses for 
the H envelope, a thick one (g(H) = 10~*) and a thin one 
(g(H) = 10"^"). We have verified that the initial masses are 
indeed sensitive to the cooling sequences used, as expected. 
We have obtained larger initial masses when considering a 
thin envelope, due to the longer cooling times obtained in 
this case. As previously pointed out, in principle it can be 
expected that the cooling time scale should be smaller for a 
thin H envelope model, b ut this assumption is only true a t 
low enough luminosities l|Prada Moroni fc Stranierol 120021 '). 
In fact, at intermediate luminosities a white dwarf with a 
thinner H envelope evolves slower than the thicker counter- 
part because it has an excess of energy to irradiate (Tassoul, 
Fontaine & Winget, 1990). The maximum difference in the 
initial masses (~ 1 Mq) has been found to occur for high- 
mass progenitors (M > 5 Mq), while this value is one order 
of magnitude smaller for smaller masses (~ 0.1 Mq). How- 
ever, it should be noted that many other combinations are 
possible, for instance, with different thicknesses of the He en- 
velope, which in this case has been kept fixed. However, since 
it is impossible to know which is the real chemical stratifi- 
cation of the outer layers of each individual white dwarf we 
have not formally introduced this error in the calculations, 
since in some cases we would be overestimating the error of 
the initial masses derived in this work. 

3.1.2 Composition of the core 

It should be taken into account that besides a CO core, white 
dwarfs can have other internal compositions. Those white 
dwarfs more massive than 1.05 Mq are thought to have a 
core made of ONe (Garcfa-Berro, Ritossa & Iben, 1997; Ri- 
tossa, Garcfa-Berro & Iben, 1996), while those with masses 
below 0.4 Mq have an He-core. ONe white dwarfs cool 
faster than CO or He white dwarfs because th e heat capacity 
of O and Ne is smaller than that of C or He (jAlthaus et al.l 
l2007h . On t he contrary. He whit e dwarfs are the ones that 
cool slower (|Serenelli et al.|[2002l '). Thus, those white dwarfs 
studied here with masses near the limits between different 
populations would be introducing an uncertainty in the cool- 
ing times obtained, since their cooling timescales are com- 
pletely different from one composition to another. For exam- 
ple, a 1 Mq ONe white dwarf cools 1 .5 times faster than a 
CO white dwarf with the same mass (|Althaus et al.l 120071 ). 
Thus, if an observed white dwarf has indeed a ONe core in- 
stead of the typical CO one, its progenitor lifetime would be 
underestimated in our analysis, and as a consequence, the 
initial mass derived would be more massive than the real 
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one. However, the exact impact of this depends also on the 
total age of the white dwarf. The smaller the total age, the 
higher the effect of considering a wrong internal composi- 
tion. 



3.1.3 Mass determinations when Tctt ^ 12000 K 

The errors reported in our study for the final masses only 
take into account the errors in the determination of the at- 
mospheric parameters, which can be derived with accuracy if 
high signal-to-noise spectra are acquired. However, it should 
be noted that this accuracy decreases considerably at low 
effective temperatures. According to Bergeron, Wesemael & 
Fontaine (1992), the atmospheres of DA stars below 12 000 
K could be enriched in He while preserving their DA spec- 
tral type. This He is thought to be brought to the surface as 
a consequence of the developement of a H convection zone. 
Depending on the efficiency of convection the star could still 
show Balmer lines, instead of being converted into a non-DA 
white dwarf. Thus, the assumption of an unrealistic chem- 
ical composition could have a large impact on the cooling 
times estimated, but mainly at low effective temperatures. 
Nevertheless, a very large fraction of the stars in our sample 
(~95 per cent) have temperatures well above this limit. 




2.0 4.0 6.0 

Mi (MJ 

Figure 2. Final masses versus initial masses for the white dwarfs 
in our sample. Solid and dotted lines correspond to weighted least- 
squares linear fits of the data. 



3.1.4 Total age of the white dwarfs 

As already pointed out, the derived initial masses depend 
on the cooling times, the total ages, the metallicity and fi- 
nally, on the stellar tracks used. Among these parameters 
the largest source of error is due to the uncertainty in the 
total ages of the white dwarfs. For white dwarfs in open clus- 
ters, the age can be usually derived with high accuracy from 
model fits to the turn-off location in a colour- magnitude dia- 
gram. The uncertainty on the age of a cluster is a systematic 
effect for stars belonging to the same cluster, since all the 
initial masses will be shifted together to larger or smaller 
masses in the final versus initial masses diagram l|WilIiamsl 
I2OO7I ). On the contrary, in the case of white dwarfs in com- 
mon proper motion pairs, the accuracy in the total age de- 
pends on the evolutionary stage of the companion. The ac- 
curacy of the age using isochrone fitting could be high if 
the star is relatively evolved and located far away from the 
ZAMS. Using the X -ray luminosity method, described in 



[Catalan et all l|2008l ). the ages derived could also be quite 
precise (from 8 to 20 per cent) if the star is relatively young 
{t s; 1 Gyr). 



3.2 The semi-empirical relationship 

Following closely recent works on this subjec t 
l|Ferrario et all l2005l : IWiUiamd l2007l : iKalirai et all |2008| ). 

we assume that the initial-final mass relationship can 
be described as a linear function. We have performed a 
weighted least-squares linear fit of the data, obtaining that 
the best solution is 



Mi = (0.117 ± 0.004)i\/i + (0.384 ± 0.011) 



(1) 



In Fig. [2] we represent all the data that we have re- 
calculated and this linear fit (dotted line). In past works, 
since there was not available data in the region of low- 
mass white dwarfs, a least-squares line ar fit led to an un- 
constrained result (jFerrario et al.ll2005l ). For this reason, a 
ficticious anchor point at low masses was used to represent 
the typical white dwarf mass of Ml ~ 0.57 Mq (jKepler et all 
[2007.y In our case, this is not necessary since we are now re- 
producing this well-established peak of the field white dwarf 
mass d istribution thanks to the new data in t he low-mass 
region l|Kalirai et al.ll2008l : [Catalan et allbOOSl ). As can be 
seen in Fig. [5] the theoretical initial-final mass relationship 
can be divided in two different linear functions, each one 
above and below 2.7 Mq, with a shallower slope for small 
masses probably due to the smaller efficiency of mass loss. 
Taking this into account we have performed a weighted least- 
squares linear fit for each region, obtaining 



Ml = (0.096 ± 0.005)Mi + (0.429 ± 0.015) 



(2) 



where the errors are the standard deviation of the coeffi- 
cients. 



for Mi < 2.7 Mq, whereas for Mi > 2.7 Mq we obtain; 

Ml = (0.137 ± 0.007)Mi + (0.318 ± 0.018) (3) 

In these expressions the errors are the standard devia- 
tion of the coefficients. These two independent fits, which are 
represented as solid lines in Fig. [2j seem to reproduce better 
the observational data than a unique linear fit (dotted line). 
Taking into account the scatter of the data and the values 
of the reduced these fits (7.1 and 4.4, respectively) 

we consider that the errors associated to the coefficients are 
underestimated. A more realistic error can be obtained com- 
puting the dispersion of the derived final masses, which is 
of 0.05 Mq and 0.12 Mq respectively. These are the errors 
that should be associated to the final mass when using the 
expressions derived here — Eqs. ([2| and Q, respectively. 
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3.3 Dependence on different parameters 

As already mentioned in the introduction, there are other 
parameters besides the mass of the progenitor that may have 
an impact on the final masses of white dwarfs (e.g. metal- 
licity or rotation). A detailed analysis of the results that 
we have obtained so far allows us to give some clues on the 
dependence of the initial-final mass relationship on these 
parameters. We discuss them below. 

3. 3. 1 Metallicity 

The sample of white dwarfs studied here covers a range of 
metallicities from Z = 0.006 to 0.040. From a theoretical 
point of view it is well established that progenitors with 
large metallicit y produce less ma s sive w hite dwarfs — see 
the relations of iDominguez et ahl (Il999l ) plotted in Fig. [1] 
Thus, one should expect to see a dependence of the semi- 
empirical data on metallicity. Our purpose in this section is 
to compare data with the same and different metallicity and 
evaluate if the differences in the derived masses are smaller 
or greater in this two cases. Open clusters are appropriate for 
carrying out such comparison, since all the stars belonging to 
a particular cluster have the same metallicity. For instance, 
in the case of the two stars from Praesepe (open circles in 
Fig.[T]) with initial masses around 3.0 Mq, the white dwarfs 
differ in AAff = 0.3 Mq. However, for the rest of stars in 
this cluster, the dispersion is significantly smaller (a factor 
of 2). Thus, the large spread could be explained by consid- 
ering that may be these objects are field stars. Ifowever, 
the sample of white dwarfs in the Praesepe cluster h as been 
weU studied (|Claver et al.ll200ll : [Pobbie et allbOO^ ) and it 
is unlikely that these stars do not belong to the cluster. 
The Hyades also contain a large number of white dwarfs, for 
which the initial and final masses have been derived with ac- 
curacy. In this case, all the data po ints fall in the reg i on lirn - 
ited by the theoretical relations of lDominguez et al.l l|l999f ). 
and the scatter is AA/f =0.1 Mq, which is the same as 
the difference between the theoretical relations correspond- 
ing to Z — 0.008 and Z = 0.02. This is the minimum scatter 
found in the observational data, since the points correspond- 
ing to the Hyades, together with those of Praesepe, are the 
ones with smaller error bars. In any case, this scatter, al- 
though smaller than in other cases, is still larger than the 
uncertainties, which prevents to derive any clear dependence 
on metallicity. As previously pointed out, the data of M37 
present the largest scatter in the sample (AMf — 0.5 Mq), 
but it should be noticed that the errors in the final masses 
are considerably larger than for the rest of the clusters. 

It is important to evaluate if the dispersion increases 
when data with different metallicities is considered. Compar- 
ing points with the same initial mass but difi^erent metallici- 
ties, it can be noted that AA/f is of the same order as in the 
case of equal metallicities. For example, comparing the data 
of M37 and the data of the Hyades — which have metallic- 
ities Z = 0.011 and Z — 0.027, respectively — for an initial 
mass around 3 A^q, AMf is ~ 0.5 Mq, but this scatter is the 
same when we compare data from M37 only. However, the 
scatter decreases when we analyse larger initial masses. If the 
data of M35 and M37 are compared (metallicites Z = 0.011) 
with the data of NGC3532 and one of the common proper 
motion pairs located at ~ 4.0 Mq (with Z — 0.02), it can 
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Figure 3. Correlation between final masses and metallicity. 



be seen that the maximum AMf is 0.2 Mq, although this 
is also what we obtain when we compare only the data of 
M35 located in this region. In the region of small masses 
we have data ranging from Z = 0.008 to Z = 0.04, but the 
maximum scatter is the same regardless of the metallicity 
of the stars (AMf = 0.1 Mq). Finally, at the high-mass 
end we find the same dispersion, although in this case the 
data have the same metallicity. Thus, considering the cur- 
rent accuracy of the available observational data we do not 
find any clear dependence of the semi-empirical initial-final 
mass relationship on metallicity. 

It is interesting to perform a more quantitative study 
of the correlation between the final masses and metallicity. 
In Fig. 13] we plot the differences between the observed fi- 
nal masses and the final masses obtained using Eqs. ^ and 
((3]) as a function of metallicity. The dashed line corresponds 
to the hypothetical case in which there is no difference be- 
tween the observational and the values predicted by these 
relations. In order to quantify the correlation in the sample 
of points presented in Fig. |3] we have calculated the Spear- 
man rank correlation coefficient. The value obtained is 0.036, 
which is very close to zero, indicating that there is an ex- 
tremely weak positive correlation between the difference in 
the final masses and metallicity. Since the Spearman corre- 
lation does not take into account the errors of the values 
considered, we have carried out a bootstrapping in order to 
evaluate the actual uncertainty on the correlation coefficient 
without any assumption on the error bars. This consists on 
choosing at random TV objects from our sample (which has 
also A*' objects) allowing repetition, and then calculating the 
new correlation coefficient for each of these new samples. We 
have performed this a large number of times (5000) obtain- 
ing a mean correlation coefficient of — 0.002±0.128. Thus, we 
conclude that the final masses and metallicity of this sample 
does not present any correlation, and that the scatter in the 
distribution in Fig. [1] is not due to the effect of metallicity. 
Of course, it should be taken into account that the observed 



The initial-final mass relationship of white dwarfs revisited 9 



final masses have been derived using the atmospheric pa- 
rameters reported by different authors who have considered 
also different white dwarf models, although the prescri ption s 
used in the fits are usually those of Bergeron et ah! (Il992l ) 
in all the cases. 

It is worth compar i ng ou r results with other works. For 
instance, iKalirai et al.l (|2005t ) claimed that they had found 
the first evidence of a metallicity dependence on the initial- 
final mass relationship. They noticed that half of their data 
of M37 (also plotted in Fi g. [T]) wer e in ag reement with the 
theoretical relationship of iMarigol l|200ll ). and considered 
this result as an indication of dependence on metallicity. 
On the other hand, in a recen t revision of the semi-empirical 
initial-final mass relationship, IWiUiamsl (|2007l ) analyzed part 
of the cluster data discussed here by deriving a binned semi- 
empirical initial-final mass relationship. This consisted in 
associating an initial and a final mass for each cluster by 
calculating the mean of the initial and final masses of the 
individual white dwarfs belonging to that cluster. Then, they 
compared these values as function of metallicity and, as in 
our case, they did not find a clear dependence of the semi- 
empirical data on metallicity. In fact, their conclusion was 
that metallicity should affect the final mass only in 0.05 Mq, 
considering an initial mass of 3 Mq , which is in good accord 
with our results. 



3.3.2 Rotation 



According to iDomi'nguez et ahl (|l996l V fast rotating stars 
produce more massive wh ite dwarfs than slow rota ting stars. 
The models calculated bv lDomfnguez et al.l l| 19961 ) including 
rotation predict that for a fast rotating star with an initial 
mass of 6.5 Mq the white dwarf produced has a mass in the 
range 1.1 to 1.4 Mq, w hich is considerably larg er than when 



^ ably larg er 

rotation is disregarded (|Dommguez et aLll999l ). Among our 



sample of white dwarfs in common proper motion pairs there 
are two stars (WD1659-531 and WD1620-391) that might 
exemplify the effect that rotation may have in stellar evolu- 
tion. The companion of W pi659-531, HD153580, is a fast 
rotating star according to iReiners fc Schmitd (|l993l '). with 
a tangential velocity v sin i = 46 ±5 km s-^ Thus, we can 
hypotetically assume that the progenitor of WD1659— 531 
was also a fast rotator. If we compare the masses derived for 
these stars (Table[l| we can note that starting from an initial 
mass of 1.58 Mq, which is more than two times smaller than 
the progenitor of WD1620-391 (3.45 Mq), it ends up as a 
white dwarf with approximately the same mass, 0.66 Mq. 
This indicates that the progenitor of WD1659— 531 lost 
less mass during the AGB phase than the progenitor of 
WD1620— 391. These differences may not be related to 
metallicity, since both progenitors had solar composition. 
Thus, we think that this could be the first evidence showing 
that rotation may have a strong impact in the evolution of 
a s tar, leading to more ma ssive white dwarfs, as suggested 
bv lDommguez et al.1 (|l996l ). 



been the subject of studies to investigate their magnetic na- 
ture. In particular, WD 0837-1-199 (also known as EG61), 
which belongs to Praesepe, is the only known magnetic 
white dwarf in an open cluster. The magneti c field of this 
star i s of 3 MG according to the study of iKawka et al.l 
(l2007f ). The mass that we have obtained for this star is 
0.81 Mq ±0.03 Mq (see Table[l]), rather large, although a 
bit smaller than the typical mass of magnetic white dwarfe, 
which is around 0.93 Mq (fWickramasinghe fc Ferrarig 
I2OO5I ) . Regarding white dwarfs in common proper motion 
pans, iKawka et al] (120071 ') also obtained circularly polarized 
spectra of WD0413-017, WD1544-377, WD1620-391 and 
WD1659— 531, finding evidence of magnetism only for the 
first of these stars. This resul t is in good agreement with 
the findings of oth er authors l|Aznar Cuadrado et al. I I2OO4I : 
[Jordan et a"Lll2007l ). In the case of WD0413-017, more com- 
monly known as 40 Eri B, the magnetic field is rather weak, 
2.3 kG, and the mass that we have derived is 0.54±0.02 Mq. 
Although the mass of this star is well below the typical 
mass of magnetic white dwarfs, it is in good agreement 
with the rest of white dwarfs studied bv lKawka et al.l (|2007l ) 
at the kG level. If we do not consider WD 0837-1-199 and 
WD0413— 017 in the fit carried out in the last section we 
obtain negligible changes on the expressions derived. 

Spectropolarimetric surveys of white dwarfs have sug- 
gested that there is a decline in the incidence of magnetism 
of stars with fields B < 10® G, although this incidence 
see ms to rise again when the field is much lower, i? < 100 
kG l|Wickramasinghe fc FerrariolbOOSh . The sample of white 
dwarfs that we have considered in this work contains two 
magnetic white dwarfs, one with a strong magnetic field and 
the other with a rather weak magnetic field. Although the 
influence of the magnetic field in stellar evolution has not 
been yet established, we flnd that the final masses obtained 
in these two cases are very different, being larger when the 
magnetic field is stronger. So, it is reasonable to think that 
the magnetic field could play a key role on the evolution 
of the progenitor star. However, with the current data on 
the magnetic white dwarfs belonging to the sample of white 
dwarfs used in this work it is not possible to favor one of the 
two main hypothesis regarding this issue: whether magnetic 
white dwarfs are more massive because the progenitors were 
also more massive (without any dependence on the magnetic 
field), or on the contrary, magnetic white dwarfs are more 
massive because the magnetic field had an infiuence during 
its evolution, favo ring the growth of th e core. 

According to iKawka et al.l (|2003l ). 16 per cent of the 
white dwarf population should be comprised by magnetic 
white dwarfs. So, among the sample of stars considered in 
this work, we could expect 9-10 white dwarfs to be mag- 
netic. Thus, spectropolarimetric observations of the current 
sample of white dwarfs used to define the semi-empirical 
initial-final mass relationship would be very useful to shed 
some light upon this subject. 



3.3.3 Magnetism 

One way to detect the presence of magnetic fields in white 
dwarfs is by performing spectropolarimetric observations. 
Some of the white dwarfs belonging to this sample have 



4 THE WHITE DWARF LUMINOSITY 
FUNCTION 

The white dwarf luminosity function is defined as the num- 
ber of white dwarfs per unit volume and per b olometric mag- 
nitude — see, for instance, llsern et al.1 (|l998l ): 
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Figure 5. Main sequenc e lifetime versus s tellar mass accord- 
ing to different au thors: iDommg uez et al] lll999il — D99 — 
iGirard i ct aT (200^ — G02 
and lWood. (,19921) — W92. 
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n(Mboi, T) = / 4){M)^Ij(T ~ tcooi - iprog)T-cooi dM (4) 

where M is the mass of the progenitor of the white dwarf, 
Tcooi = d,t/dM]joi is its characteristic cooling time, Mi and 
Ms are the minimum and maximum mass of the progenitor 
star able to produce a white dwarf with a bolometric mag- 
nitude Mboi at time T, tcooi is the time necessary to cool 
down to bolo metric magnitude M boi — for which we adopt 
the results of lSalaris et all (|2000l ') and lAlthaus et all (|2007l ') 
for CO and ONe white dwarfs, respectively — tprog is the 
lifetime of the progenitor, T is the age of the population un- 
der study, tp{t) is the star formation rate — which we assume 
to be constant — and 0(M) is th e initial mass fu nction, for 
which we adopt the expression of lSaloeteij (|l955l l. 

4.1 The influence of the progenitors 

To compute the white dwarf luminosity function it is also 
necessary to provide a relationship between the mass of the 
progenitor and the mass of the resulting white dwarf, that 
is, the initial-final mass relationship. Additionally, the in- 
fluence of the progenitors in Eq. ^ appears through the 
age assigned to the progenitor, which determines the star 
formation rate, and through the cooling time, fcooi and the 
characteristic cooling time, TcooI, which depend on the mass 
of the white dwarf. In order to evaluate the influence of these 
inputs, we have computed a series of theoretical white dwarf 
luminosity functions using several initial-final mass relation- 
ships (Fig. m and evolutive tracks for the progenitor stars 
(Fig. [5]). Whenever possible we have adopted the same set 
of stellar evolutionary inputs, that is, the initial-final mass 
relationship and the main-sequence lifetime corresponding 
to the same set of calculations. 



We compare the resulting theoretical luminosity func- 
tions with the data obtained by averaging the different ob- 
servational determinations of the white dwarf luminosity 
function (Knox, Hawkins & Hambly, 1999; Leggett, Ruiz & 
Bergeron, 1998; Oswalt et al. 1996; Liebert, Dahn & Monet, 
1988). The theoretical white dwarf luminosity functions were 
also normalized to the observational value with the small- 
est error bars in number density of white dwarfs, that is 
log(iV) = -3.610, log(L/L0) = -2.759, avoiding in this 
way the region in which the cooling is dominated by neutri- 
nos (at large luminosities) and the region in which crystal- 
lization is the dominant physical process — at luminosities 
between log(L/LQ) ~ —3 and —4. 

Fig. [S] shows the resulting white dwarf luminosity func- 
tions when different stellar evolutionary inputs are used. At 
low luminosities it can be noticed the characteristic sharp 
down-turn in the density of white dwarfs. This cut-off in 
the number counts has been interpreted by diff erent authors 
IWinget et al.lll987l : iGarci'a-Berro et al.lll988l ) as the conse- 
quence of the finite age of the Galactic disc. Thus, a compar- 
ison between the theoretical luminosity functions and the 
observational data can provide information about the age 
of the Galactic disc. In this figure the cut-off of the ob- 
servational white dwarf luminosity function has been fitted 
using an age of the disc of T = 11 Gyr for all the ca s es ex - 
cept for the case in which the expressions of I Wood! (|1992D 
were used. In this last case the best-fitting is obtained us- 
ing T = 10.5 Gyr. This can be understood by comparing 
the different stellar evolutionary inputs considered in this 
work. As c an be s een in Fig. U the initial-final mass rela- 
tionship of lWoodI (I1992I ) is the one that produces less mas- 
sive white dwarfs. The semi-empirical relationship that we 
have derived in this wo r k is s imilar to that of Hurley, Pols 
& Tout (2000). iMarigol | 200 ll ) predicts more massive white 
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d warfs at the intermediate mass domain, whereas the results 
of iDommguez et al.l (|l999l ) produce more massive remnants 
for the low-mass end and less massive white dwarfs for the 
high mass end, but always substantially larger than those 
obtain ed with the initial-final mass relationship of IWoodI 
These differences mainly arise from the procedure 
used to calculate each theoreti cal initial-final mass rela tion- 
ship. The rela tions derived bv lDommguez et"al] (|l999l ') and 
iMarigol { 200 ll ) are obtained using fully evolutive models but 
using different treatments of convective boundaries, mix- 
ing and mass-loss rate s on the AGB phas e. On the other 
hand, the relations of iHurlev et al.1 l|2000l) were obtained 
from a fitting of the observational data from eclipsing bina- 
ries and open clusters, obtaining analytic formu l ae for differ- 
ent metallicities. Finally, the relation of IWoodI (|l992l ) is an 
exponential expression derived from a fit to the PNN mass 
distribution. In Fig. [S] we show the different main-sequence 
lifetimes as a function of the main-sequence mass. In those 
cases in which there was a dependence on metallicity we 
have adopted Z = Zq. As it can be seen there, the be- 
haviour of the different mai n-sequence lif etimes is very sim- 
ilar, although for the case o f lWoodI (|1992| ). stars spend more 
time in the main-sequence than in the rest of cases, espe- 
cially for those stars with large masses. Considering this, one 
should expect that the fit of the cut-off of the white dwarf 
luminosity function would correspond to a longer Galac- 
tic disc age when using the expressions of IWoodI l|l992h . 
On the contrary, we have obtained a youn ger Galactic disc. 
The longer progenitors lifetimes of IWoodI ((1992) are in part 
compensated by the fact that its corresponding initial-final 
mass relationship favors the production of low-mass white 
dwarfs, which cool down faster at high luminosities. A simple 
test can be done by using, for example, the stellar tracks of 



iDommguez et ahl l|l999l ), which give shorter progeni t or life - 
times, and the initial-final mass relationship of lWoodI (|l992h . 
In this case we obtain a Galactic disc age of 10 Gyr, 1 Gyr 
younger than if we consider any of the other initial-final 
mass relations hips shown in Fig. |4] Thus, the behaviour of 
the relation of IWoodI (|l992l ) is clearly different from others 
in the literature, and this has important implications on the 
resulting white dwarf luminosity functions. 

Comparing the results of Fig. [6] it can be noticed that 
for the hot end of the white dwarf luminosity function there 
are not differences whatsoever. In fact, all the theoretical lu- 
minosity functions are remarkably coincident. The only vis- 
ible differences, although not very relevant, occur just after 
the crystallization phase has started, at log(L/LQ) ~ —4.0. 
Two essential physical processes are associated with crystal- 
lization, namely, a release of latent heat and a modification 
of the chemical concentrations in the solid phase. Both pro- 
vide extra energy sources and lengthen the cooling time of 
the star. This is the reason why all the theoretical calcula- 
tions predict a larger number of white dwarfs for these lumi- 
nosity bins. Beyond the cut-off, we find that the density of 
white dwarfs is smaller in the case of Wood (1992^), and this 
is because this region is dominated by massive white dwarfs, 
since the progenitors of these stars were also massive, and as 
shown in Fig. [S] spent less time at the main-sequence. Thus, 
massive white dwarfs have had enough time to cool down 
to such low luminosities. In any case, it is worth mentioning 
that low-mass white dwarfs cool faster at high luminosities, 
but for small luminosities it occurs just the opposite, since 
massive white dwarfs crystallize at larger luminosities, and 
this implies smaller time delays. 



4.2 The luminosity function of massive white 
dwarfs 

4-2.1 Effect of the initial-final mass relationship 

The influence of the initial-final mass relationship on 
the white dwarf luminosity function should be more ev- 
ident when it is con strained to massive white dwarfs 
(|Di'az-Pinto et al.lll99j '). Recently, Liebert, Bergeron & Hol- 
berg (2005b) performed high signal-to-noise spectroscopic 
observations of more than 300 white dwarfs belonging to 
the Palomar Green (PC) Survey. The analysis of this set of 
data has provided us with a sample of white dwarfs with 
well determined masses that allows for the first time the 
study of the white dwarf lumi nosity function of massive 
wh ite dwarfs (|lsern et al] |2007| ) . Unfortunately, according 
to iLiebert et al.1 (|2005bl ) the completeness of the sample 
decreases severely near 10 000 K, where white dwarfs with 
small masses (0.4 Mq) are brighter (Mv ~ 11) than massive 
white dwarfs with M > 0.8 Mq (which have visual magni- 
tudes around Mv ~ 13). Consequently, above My = 11 this 
survey only has detected white dwarfs with masses larger 
than 0.4 Mq. Therefore, we will limit the analysis to white 
dwarfs brighter than Mv ~ 11. 

We have computed a set of white dwarf luminosity func- 
tions considering an age of 11 Gyr for the Galactic disc 
and using bins of visual magnitude. In Fig. [7] we show from 
top to bottom the total luminosity function and the lumi- 
nosity functions of white dwarfs with masses larger than 
0.7 Mq and 1.0 Mq, respectively. The total luminosity func- 
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Figure 7. White dwarf luminosity functions versus vi- 
sual magnitude using d ifferent initial-final mass relationships: 
Domi'ngu ez et alj l ll999l) — D99 — iMarigol l l200lt) — MOl — 
Wood (199^ — W92 — and this work. From top to bottom we 
show the total luminosity function, and the luminosity functions 
of white dwarfs with masses larger than 0.7 Mq and 1.0 A/q. Cir- 
cl es, triangles and squa res correspond to the observational data 
of lLiebert et al.l ll2005bh . 
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Figure 8. White dwarf luminosity functions versus visual magni- 
tude using the initial-final mass relationship derived in this work 
and different star formation rates. Circles, triangles and squares 
correspond to the observational data of iLiebert et al.. (.2005b.) . 



Figure 9. Same as Fig. [8] but using the initial-final mass rela- 
tionship of Wood (1992). 



tion (that is, considering the whole the range of masses) 
was normalized to the bin corresponding to A/y = 11, and 
then, this normalization factor was used for the luminosity 
functions of white dwarfs more massive than 0.7 Mq and 
1.0 Mq. In this case, we have us ed the stellar evolution- 
ary iri puts oflP^ommeuc z et all |l999). Marigo (2001), Wool 
(il992h and the semi-empirical initial-final mass relationship 
that we have derived in the previous section. Comparing the 
different theoretical luminosity functions, it can be noted 
that the predicted number of massive white dwarfs i s larger 
when using the inputs of iDommguez et ahl (|l999l ) , iMarigol 
(I2OO1I ) and our semi-empirical initial-final mass relationship 
in comparison wi th the result s obtained when considering 
the expressions of lWoodl(ll992l ). This is obviou sly due to the 
fact that the initial-final mass relationship of 'Wood' (199^ 
favors the production of low-mass white dwarfs. It can also 
be noted that the density of massive white dwarfs is slightly 
la rger wh e n con sidering the initial-final mass relationship 
of lMarigol (1200 ll ) than that obtained when using the initial- 
final mass relationship derived here. The reverse is tru e when 
the initial-final mass relationship of lDomfnguez et al.l (|l999l ) 
is used. With out considerin g the results obtained when the 
expressions of lWoodI l|l992l ) are used, it can be noted that it 
is not possible to evaluate which initial-final mass relation- 
ship produces a theoretical luminosity function that better 
fits the observational data, since the error bars of the ob- 
servational data are larger than the differences between the 
theoretical results. In any case, what it can be clearly seen is 
that all the theoretical relations predict more massive white 
dwarfs than the observations when a mass cut of 1.0 Mq is 
adopted, except in the case of lWoodI l| 19921 ). 



^.2.2 Effect o] the star formation rate 



The star formation rate considered in our calculations has 
also an important influence on the number of massive white 



The initial- final mass relationship of white dwarfs revisited 13 



dwarfs produced. In our previous calculations we have con- 
sidered a constant star formation rate, which heis led us to 
obtain more massive white dwarfs than those detected in the 
PG Survey. To evaluate the effect of the star formation rate 
on the white dwarf luminosity function we have repeated the 
calculations considering the semi-empirical initial-final mass 
relationship derived in this work and an exponentially de- 
creasing star formation rate tp{t) = exp(— f/r), with r = 3 
and 5 Gyr, respectively. As it can be noted from Fig. |8]the 
production of massive stars decreases considerably when a 
exponentially decreasing star formation rate is considered. 
The observational data corresponding to the whole range of 
white dwarf masses and to masses larger than 0.7 Mq is 
better fitted when a constant star formation rate (r oo) 
is assumed in the theoretical calculations. On the contrary, 
for white dwarfs with masses larger than 1.0 Mq the agree- 
ment is better if we consider a variable star formation rate. 
However, for such massive stars there is only data for one 
magnitude bin. More observations corresponding to massive 
white dwarfs are needed to confirm this behaviour of the 
luminosity function. 

For the sake of comparison we have carrie d out the sam e 
calculations but considering the expressions of lWoodl(|l992h . 
The resulting luminosity functions are shown in Fig. [5] In 
this case, the observational data for stars with masses larger 
than 0.7 Mq is not fitted regardless of the star formation 
rate assumed. On the contrary, the fit is better when consid- 
ering stars with masses above 1.0 Mq. Since there are more 
observational data for the range of masses above 0.7 Mq, 
we find more reliable the conclusions that can be obtained 
from a comparison of these data than those obtained from 
only one magnitude bin, as is the case of more massive white 
dwarfs. Thus, it seems clear th at obse r vation s rule out the 
initial-final mass relationship of lWoodI l|l992l ). and that the 
other initial-final mass relationships used in this work seem 
to be more reliable, although the present status of the ob- 
servational data does not allow to draw a definite conclusion 
about which initial-final mass relationship is more adequate. 



5 THE WHITE DWARF MASS 
DISTRIBUTION 

The understanding of the precise shape of the mass distri- 
bution of white dwarfs off'ers a sorely needed insight about 
the total amount of mass lost during the course of the last 
phases of stellar evolution. Thus, a detailed study of this 
function, from both the theoretical and the observational 
perspective s can give us clues on the initial-final mass re- 
lationship l|Ferrario et al.|[2005h . For that purpose, we have 
computed a series of theoretical mass distributions using 
different evolutive stellar models and their corresponding 
initial-final mass relationships, if available. As in the case of 
the white dwarf luminosity function we have adopted an age 
of 11 Gyr for the Galactic disc, a constant star formation 
rate and the initial mass function of Salpeter. The theoreti- 
cal mass distributions were then normalized to unit area. 

Our purpose is to c ompare our results wit h the observa- 
tional data obtained bv lLiebert et al.l (l2005bl ) from the PG 
Survey — left panel of Fig. 1101 — and the recent data ob- 
tained by iDeGennaro et~aL I (|2008D from the SDSS — right 
panel of Fig.[lOl As previously pointed out, the accuracy on 



the mass determinations decreases considerable when white 
dwarfs are cooler than 12 000 K. Hence, both the theoreti- 
cal and the observational mass distributions in this section 
consider white dwarfs with T^a ^ 12 000. In the case of 
the data from the PG Survey, we have computed the ob- 
servational mass distribution con sidering the Vlnax values 
reported by (|Liebert et al.l [20051J ) taking into account the 
errors in the masses assuming a gaussian distribution. Then, 
the final distribution has been normalized to unit area, as 
it was done when considering the theoretical mass distribu- 
tions. In the case of the data from the SDSS, the observa- 
tional data s hown in Fig. 1 10 | (right) is the mass distribution 
computed bv IDeGennaro et al.l (|2008h . 

In Fig.[Tn]we pl ot our results cons i dering t he stellar evo- 
lutionary i nputs o f iDommguez et "all (|l999l ), iGirardi et al.l 
(2002) and IWoodI (| 19921 ). We also show our results when 
considering the evolutive stellar models of iDommguez et al.l 
(Il999|) and the initial-final mass relation derived in 



this work, a s well as the relation recently obtained by 
iKalirai et al.l l|2008l ). All the white dwarf distributions have 
been normalized to the total density obtained in each case. 
As it can be noted, there is a well defined peak in all the 
mass distributions, the location of which is defined mainly 
by the initial-final mass relationship considered. On the con- 
trary, the height of the peak depends also on the lifetime 
of the progenitors. This can be understood with the help of 
Fig. [l] The most abundant stars a r e thos e wit h small masses 
2^1 Mq). In the case of I Wood! l|l992l) and iGirardi et all 
2002h . the white dwarfs corresponding to these progenitors 
have masses well below ~ 0.6 Mq, and this is the reason why 
the central peak is located at smaller masses in Fig. [10] This 
peak is then shifted to larger masses when the initial-final 
mass relationship considered favors the production of more 
massive white dwarfs for the low mass progenitors. This is 
th e case of the s e mi-em pirical initial-final mass relationship 
of iKalirai et al. ll2008h. our relatio nship or the theoretical 
relation of Dommguez et all (|l999l ) . Note that in the latter 
case, the production of massive white dwarfs is not favored, 
opposite to what occurs to our semi-empirical relationship, 
but the peak is located at the same mass since low-mass 
progenitors are the ones that dominate. 

If we compare th e theoretical mass fun ctions with the 
observational data of iLiebert et al. I l|2005bD from the PG 
Su rvey (left panel of Fig. | 1 0l) and the recent data obtained 
bv I DeGennaro et al.l \2Qm from the SDSS (right panel of 
Fig. I10[> it can be noted how in both cases the location of the 
central peak is well fitted by the predictions corresponding to 
o ur semi-empirical relation ship and the theoretical relation 
of lDonifnguez et al.l (|l999 ^. The height of the peak is better 
fitted when considering our semi-empirical relationship, al- 
though a bit lower in comparison with the SDSS data. The 
sample of white dwarfs corresponding to the SDSS is very 
complete since it includes the 1733 stars with <; ^ 19 and 
Teft ^ 12 000 K present in the SDSS DR4, which is seven 
times the population covered by the PG Survey. 

It is worth mentioning that we have not included the 
He white dwarf population in our calculations. This is the 
reason why there are no white dwarfs below a certain mass 
threshold ( 0.45 Mq in the case of our semi-empirical initial- 
final mass relationship). The observational data (both from 
the PG survey or the SDSS) do indeed present a certain 
number of white dwarfs in this low-mass region. On the 
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data collected in the PG Survey (left) and those obtained bv iDeGennaro et al., (,2008i ) corresponding to the data in the SDSS (right). 



other hand, we have included the ONe white dwarf popula- 
tion which is thought to be dominant for masses larger than 
1.05 Mq. ONe-core white dwarfs cool faster than CO-core 
white dwarfs due to the large heat capacity of C in compar- 
ison with that of O and Ne. For this reason a bump located 
at 1.05 Mq in the mass distribution that corresponds to 
the change of the cooling rate is clearly visible. Since the 
cooling rate is larger for the ONe white dwarf population, 
there is an increase of the number of white dwarfs produced 
in this region. As it can be noted in Fig. [10] the density of 
massive white dwarfs is remarkably l ower w h en co nsidering 
the initial-final mass relationship of IWoodI l|l992l ) since it 
predicts less massive white dwarfs. In the rest of cases the 
white dwarf mass distributions are approximately coincident 
in this region. 



6 SUMMARY AND CONCLUSIONS 

In this work we have revisited the initial-final mass relation- 
ship and discussed it in a comprehensive manner. With this 
purpose, we have re-evaluated the available data in the lit- 
erature, mainly based on open clusters, that are currently 
being used to define the semi-empirical initial-final mass re- 
lationship. We have used the atmospheric parameters, total 
ages and metallicities report ed in the litera t ure an d followed 
the procedure described in I Catalan" et all (|2008l ) to derive 
the initial and final masses of these white dwarfs. Thanks 
to these data and our own work based on common proper 
motion pairs we have been able to collect a very hetero- 
geneous sample of white dwarfs, covering a wide range of 
ages and metallicites. Most importantly, with this study we 
have covered the range of initial masses from 1 to 6.5 Mf^ , 
which was poorly covered until recently (|Kalirai et al.|[200a : 



[Catalan et al.|[2008l ). The extension of the initial-final mass 
relationship to the low-mass end is important since we are 
now studying the most populated region of initial masses 
according to the initial mass function of Salpeter, and this 
means that we are also reproducing the we ll-established 
peak of the field white dwarf mass distribution er et al.l 
l2007f ). 

As discussed previously, for each cluster the results 
present an intrinsic mass spread, which may indicate that 
mass loss could depend more on individual stellar proper- 
ties than on a global mechanism. Thus, mass-loss processes 
could even be a stochastic phenomenon (Reid 1996), being 
thus impossible to reproduce with models. Another expla- 
nation for the spread of masses found for each cluster could 
lie on the fact that we do not know the internal composition 
of the white dwarfs in our sample, and we are using cooling 
sequences that have fixed values for the C/O ratio at the 
core and most importantly, fixed thicknesses of the H and 
He envelopes which may not be appropriate to describe the 
cooling of individual white dwarfs in some cases. This fact 
could affect the cooling times d erived in approximately 10 
per c ent at very low luminosities (jPrada Moroni fc Stranierd 
|2002|) . 

Since there is no compelling reason to justify the use 
of a more sophisticated relationship, we have performed a 
weighted least-squares linear fit of these data, and from a 
detailed analysis we have found that there is no correla- 
tion in this sample between the final masses and metallicity. 
We have also given some clues on the dependence of the 
initial-final mass relationship on other parameters, such as 
rotation and magnetism. In the case of rotation, we have 
probably found the first evidence that rotation may have 
in stellar evolution, corroborating that when the progeni- 
tor is a fast rotating star, the resulting white dwarf is more 
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massiv e, in agreement with the study of iDommguez et alj 
1 19961 ). Among our sample there are two magnetic white 
dwarfs with rather different magnetic fields. The masses de- 
rived in both cases are also clearly different, being more mas- 
sive the one with stronger magnetic field. This indicates that 
the intensity of the magnetic field might be related to the 
mass of the white dwarf produced. However, given the lack 
of information on other possible magnetic white dwarfs be- 
longing to this sample, our results are not conclusive enough 
to assess the impact that magnetic fields may have in the 
structure of white dwarfs. 

In the second part of this work we have tested the 
initial-final mass relationship by studying its effect on the 
luminosity function and mass distribution of white dwarfs. 
For this purpose we have used different stellar evolutionary 
inputs (stellar tracks and initial-final mass relationships). 
We have also computed the luminosity function of massive 
white dwarfs in order to evaluate the impact of the initial- 
final mass relationship. We have noted some differences be- 
tween the theoretical luminosity functions, obtaining a clear 
dependence on the considered initial-final mass relationship, 
as expected. From a comparison of our results with the ob- 
servational data from the Palomar Green Survey we have 
always obtained a reasonable fit when the range of masses 
was constrained to M > 0.7 M0, ex c ept w hen using the 
initial-final mass relationship of lWoodI l| 19921 ). These calcu- 
lations were performed assuming a constant star formation 
rate, but we have also computed the luminosity functions 
considering an exponentially decreasing star formation rate. 
We have shown that the production of massive white dwarfs 
is dependent on the assumed star formation rate, since in the 
latter case the density of massive white dwarfs drops con- 
siderably. Given the presently available observational data, 
any attempt to discern which initial-final mass relationship 
better fits the data is not feasible. However, our results favor 
an initial-final mass relationship th at produces m ore massive 
white dwarfs than the relation of IWoodI l|l992l ). This is an 
i mportant find ing, since the initial-final mass relationship 
of lWoodI (| 19921 ) is the most commonly used relationship for 
computing the theoretical white dwarf luminosity function. 

In the case of the white dwarf mass distribution we 
have obtained more conclusive results. From a compari- 
son of our results with th e observational data available 
from both the PG survey jLiebert et all l2005bl ) and the 
SDSS jPeGennaro et al.ll2008l ) we have noted that the semi- 
empirical initial-final mass relationship derived in this work 
is the one that better fits the central peak of the mass dis- 
tribution, being in this way very representative of the white 
dwarf population. On the contrary, the agreement w ith the 
observ atio nal data is rath e r poo r when the relations o f lWoodI 
1 19921 ) and lGirardi et all (|2002l ) are considered. 

The study carried out in this work evidences the neces- 
sity of increasing the number of high quality observations of 
white dwarfs belonging to stellar clusters, or common proper 
motion pairs, or to any system that may allow the determi- 
nation of their total ages and original metallicities with ac- 
curacy. On the other hand, this increase in the observational 
data should be accompanied by refined theoretical studies 
of the initial-final mass relationship. 
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